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We use the "modes of the universe" approach to study a cavity-mediated two-photon logical 
gate recently proposed by Koshino, Ishizaka and Nakamura. We clarify the relationship between 
the more commonly used input-output formalism, and that of Koshino et al., and show that some 
elements of this gate had been anticipated by other authors. We conclude that their proposed gate 
can work both in the good and bad cavity limits, provided only that the pulses are long enough. 
Our formalism allows us to estimate analytically the size of the various error terms, and to follow 
the spectral evolution of the field + cavity system in the course of the interaction. 

PACS numbers: 03.67.Lx, 42.50.Ex, 42.50.Dv 



I. INTRODUCTION 



In a recent publication K. Koshino, S. Ishizaka and 
Y. Nakamura introduced a cavity-mediated scheme to 
implement a deterministic photon-photon VSWAP gate 
which can be operated, in principle, in a completely pas- 
sive way: no external pulses or fields are required either 
to initialize or read-out the atom (or equivalent three- 
level system) in the cavity, nor to manipulate its inter- 
nal state or energy level structure in between the single- 
photon pulses. This could represent a substantial simpli- 
fication over existing cavity-based proposals for single- 
photon quantum logic, which have followed on the pio- 
neering work Q of L. M. Duan and H. J. Kimble. 

Our goal in this paper is to fully characterize the condi- 
tions under which the proposal by Koshino et al. works, 
and to clarify its relationship to other previously-known 
results. We choose do this by using the "modes of the uni- 
verse" formalism [!, H[ , in which the "quasimodes" of the 
optical cavity are written as superpositions of modes of a 
much larger cavity (the "universe") that encloses it. This 
formalism is closely related to the one used by Koshino 
et al., and so we can use it to show that their results are 
in fact consistent with expressions obtained for related 
systems, over the years, by other workers making use of 
the simpler Collet- Gardiner input-output formalism 
after accounting for a non-obvious phase factor difference 
between the two approaches. We deri ve expre ssions for 
the detuning needed to perform the VSWAP gate that 
arc more general than the ones in [l| , being also valid in 
the good-cavity limit (g > k, where g is the atom-cavity 
field coupling, and k the bare cavity decay rate). 

Our formalism allows us also to obtain simple analyt- 
ical estimates for the "error terms" arising from the fail- 




FIG. 1: The setup considered in this paper. The small 
(I <C L) one-sided cavity has a partly reflecting mirror with 
amplitude reflection and transmission coefficients r and t 
(f 2 + P = 1), and contains a single three- level system with 
two degenerate transitions, corresponding (for example) to 
horizontal and vertical polarizations H and V . 



ure of various approximations, or from competing pro- 
cesses such as spontaneous emission, and provides a nat- 
ural framework for numerical calculations. We present 
results of these calculations, including the evolution of 
the field + cavity system in the course of the interaction. 



II. THE "MODES OF THE UNIVERSE" 
FORMALISM 

We use the formalism developed in the paper Q (es- 
pecially Appendix A), which describes a leaky cavity 
bounded by a perfect mirror at z — I and a semitranspar- 
ent mirror at z = 0, and an auxiliary cavity (the "out- 
side world") bounded by a perfect mirror at z = — L 
(L — > —oo) and the mirror at z — (see Figure 1). 

The mode functions are (Eq. (2.2) of Q) 



U k (z) 



|£fc smk(z + L) 
I M fe sinfc(z - I) 



for z < 0, 
for z > 0. 



(1) 
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The £fc are taken to be alternately +1 and —1, and the 
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Mk then are given by 



M k = 



(ck/Z) 



1/2 



[(fi fe -r> c )2 +K 2]l/2 



(2) 



Here we have changed the notation of 3 slightly, so that 
k, rather than T, is the cavity amplitude decay rate. tt c 
is the central frequency of the cavity quasimode under 
consideration. The allowed frequencies flk = ck and 
wavevectors k for the modes of the universe can be de- 
rived from the eigenvalue equation (Eq. (A2) of Q) 



fsm[k{L-l)] = sin[k(L + l)] 



(3) 



where f is the amplitude reflection coefficient of the 
cavity input mirror. The incoming field has the form 
(Eq. (2.6a) of d) 



E : 



( + ) 



On ^ 



2i \e AL 



1/2 



ikL—iQkt 



(4) 



and the cavity field is then (Eq. (2.9) of 



cav 2i^\Ae Q Al) 

k x 



1/2 



l(t)t 



with the cavity "single quasimode" operator 



a(t) 



E M k a k< 



iki-i(n k -n c )t 



(5) 



(6) 



In this work we actually need to consider two cavity 
quasimodes, corresponding to horizontal and vertical po- 
larizations. We will denote the corresponding "modes of 
the universe" operators by dhk and a V f~, and the quasi- 
mode operators by ah and a v . We describe the coupling 
of the atom to the field by the standard interaction- 
picture Hamiltonian hg(a<j^ + a^cr), where a and are 
the usual atomic lowering and raising operators. The 
Hamiltonian then takes the form 



' i5kt ~ al k aoie i5kt ) (7) 



with 



M' k = 



211/2 



(8) 



and 5k = (fife — Q c ) + (^c — ^o) = ilfe — ujq = ck — luq. 
Here, both the resonant atomic frequency ujq and the 
coupling constant g are taken to have the same values 
for the H and V transitions. This Hamiltonian is to 
be compared with Eq. (f) of The main difference 
is that we have allowed, through the coefficients ML for 
different coupling strength of the cavity to outside modes 



with different frequencies, something which is important 
in the good cavity (small n) limit. 

We note that a factor e lkl in ([6]) has been neglected 
in writing Eqs. (|7|) and ([5]). The assumption is that the 
small cavity is so small that the factor is essentially con- 
stant as a function of k, that is, if k = k c + Ak, then 
k c l ~ 2n7r (resonance condition) and Akl ~ kI/c <^ 1 
(since, as we shall see later, the range of modes that 
eventually develop appreciable amplitudes is of the order 
of k in frequency space). This assumption may only have 
to be reevaluated for a very bad cavity; note that in our 
formalism, k = ct 2 /4Z, so kI/c ~ t 2 /4, which under most 
circumstances should indeed be very small. 

Always under the assumption that I <C L, we can take 
k = k c +rvK I L, where n is a (positive or negative) integer, 
and k c — fl c /c. We do not expect to have more than one 
photon in the system at any given time, so the state of 
the system can be written as 



|V(*)> =C 2 (t)|2)|vac) ft |vac)„+^C hn (t)|l)|l) 

n 

+ £a„(*)|0}|l>Jvac) ft 



vac) 



(9) 



where the first ket refers to the state of the atom, and 
the second and third to the field; | vac) ;,(.„) is the vacuum 
state for the horizontal (vertical) modes; \l)hn is a state 
with one photon in horizontal mode n; and \l) vn is a 
state with one photon in vertical mode n. 

Using ([7]) and ([9]), one gets the equations of motion 



c 2 


n 


Chn 


= -gM' n C 2 e iSnt 




= -gM' n C 2 e^ 1 


with 





m e- iS ^ (10a) 



ML = 



(ck/L) 



1/2 



[(nnc/L) 2 + n 2 } 1 / 2 



(10b) 
(10c) 



(11) 



and 5 n = nnc/L + 5 a , where S a = Q c — wo is the atom- 
cavity detuning. 

These equations have to be integrated with initial con- 
ditions chosen so that the incoming field ([JJ represents a 
suitable pulse. As an example, the set of coefficients be- 
low, which we have used for our numerical calculations, 
describes an incoming Gaussian pulse initially centered 
at zq (with — L < zq < 0), with a duration T, that 
is to say, a spatial width cT, and a carrier frequency 
fl c + aiQir/L. The possibility of having the initial set of 
coefficients peak at a frequency different from the cavity 
frequency is the way to account for a possible detuning 
between the external field and the cavity in our formal- 
ism; here this detuning A = O — fl c — cn n/L. (The 
analytical calculations in the next Section show how un- 
der some conditions one can "pull out" A, to deal with 
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it in more conventional ways. 



«•<»)- (try?-- 



innza/L ^ — (cT(n — no)7r/2L) 



(12) 



Note that the coefficients in (fT2|) are (approximately) 
normalized to unity. That is, Sn|C™(0)| 2 = 1j P r0 " 
vided that the sum extends over a sufficiently large num- 
ber of modes. For example, in our numerical work we 
take C vn {0) = C n (0), along with C 2 (0) = and all 
Chn(0) = 0. Of course, one could alternately take 
C hn {0) = C n (0), with C 2 (0) = and all C vn (0) = 0. 

In principle, a very long initial pulse can be accurately 
represented by relatively few modes, of the order of a few 
times L/cT; however, it is important to note that in the 
course of the interaction many more modes, of the order 
of kL/c, may become appreciably excited, as we shall 
show below. 

The formalism introduced here provides a simple 
framework for numerical calculations, and we shall show 
the results of some such calculations in Section IV be- 
low, but it is worth noting that substantial progress can 
be made by analytical methods as well. Accordingly, the 
following section presents a formal solution of Eqs. (|T0|) 
that is quite general (with no further approximations), as 
well as the simpler results of an adiabatic approximation, 
valid for sufficiently long pulses. 



III. ANALYTICAL RESULTS 

A. General solution 

Clearly, Eqs. (TTOj) can be simplified by introducing two 
sets of variables C n ± = {Chn ± C vn )/y/2. Then the C„_ 
are constant, and the system ([T0|) reduces to 



C 2 =gV2j2KCn+t 

n 

C n+ = -gsplM'C^ 1 



-i8„t 



(13a) 
(13b) 



That is, the degenerate three level system is formally 
equivalent to a two-level system, coupled to the polariza- 
tion H + V, plus a "dark state" (|0) - \l))/y/2. Hence, 
we expect our results, in an appropriate limit, to reduce 
to many results already available for a single two-level 
system in a cavity, interacting with a single photon pulse 

m- 

We can formally integrate equation f|13b[) and substi- 
tute in ([T5a)l . We get 

C 2 = - 2g 2 V M' n 2 f e-W-^dtf) dt' 

+ 9 V2j2M^C n+ (0)e- u " t (14) 

n 

The second term in (TT4")) is a driving term that depends 
on the initial condition (that is, the initial shape of the 



wavepacket); below we shall call it fo(t) for brevity. In 
the first term, we can use 



AI, 



,2 C K 



and 



L (m:c/L) 2 + k' 2 



r crni 

On = —j- + "a 



(15) 



(16) 



and convert the sum over n into an integral over w, with 
uj = cmr/L: 



E 



L ^— ' (rmc/L) + K 

ti — — nri v ' ' 



^ e -i(cnir/L)(t— t') 



,-it\t-t'\ 



(17) 



So now (3) becomes (since t' < t, by construction) 



C 2 ~ -2g 2 / e-( K +^)(*- t ')C 2 (t') dt' + f (t) (18) 
Jo 



with 



f(t)=gV2'£M n C n+ (0) t 



-iS n t 



(19) 



This integro-differential equation is of a form that de- 
scribes a damped, driven harmonic oscillator and so can 
be solved exactly. To make it slightly more general, we 
can even add a term — 7C2 to the right-hand side of 
(|18[) , to account approximately for spontaneous emission 
losses. Such an approach breaks the unitarity of the sys- 
tem (which means the total probability is no longer con- 
served), and does not properly account for the fact that 
the atom must return to one of the ground states after 
a spontaneous emission event, but it should be a reason- 
able lowest-order approximation when the spontaneous 
emission probability is small. 

We proceed by taking the derivative of (|18p . with the 
optional term — 7C2, and using flTBl itself to eliminate 
the integral, with the result 

C 2 + (7 + k + i5 a )C 2 + 2g 2 + 7 (k + iS a ) - F(t) (20) 



with 



F(t) = f + ( K + id a )f 



(21) 



With the initial condition C 2 (0) = C 2 (0) = (atom un- 
excited before the pulse arrives), the formal solution of 
this equation is 



C 2 (t) = 



1 



Ai — A 2 Jo 



F(t') dt' 



(22) 
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in terms of the eigenvalues, Ai and A2, of the character- 
istic equation 

A 2 + (7 + k + iS a )X + 2g 2 + 7 (« + iS a ) = (23) 

We stress that this represents, in principle, a full solution 
to the problem, valid for any type of cavity (good, bad 
or in between) and any kind of pulse (fast or slow). One 
only needs to use Eq. ()19[) (preferably in integral form) 
to calculate f(t), then obtain C2 from Eqs. (|21||2"3"1) . then 
substitute that in Eqs. (| 13b[) . from which one can get all 
the coefficients of the state Q, and answer any question 
about the quantum logical operation performed, includ- 
ing its fidelity. Further simplification of the results is 
possible, however, in a particularly important case, that 
of a very slow (adiabatic) pulse, which we consider in the 
next subsection. 



B. Adiabatic approximation for slow pulses 

In this section we consider the case in which the time 
scale over which the incident pulse changes, given by T, 
is much longer than all the other scales in the problem, 
and, most importantly, those set by the solutions Aj. 2 of 



T 



< |Al,: 



(24) 



In this case, things can be simplified, as follows. In the 
expression (p~9|) . the frequencies 8 n will be of the form 
8 n = c(n — no)ir / L + A + 8 a , where c(n — ti^-k/L is 
expected to be a slow frequency, and we allow for the 
possibility of either or both detunings being large (recall 
A = uqCtt/L). We then write 



f(t) 



-i(A+S a )t 



fo(t) 



(25) 



where /o is a slowly-varying function, and approximate 
(1211 by 

F(t) ~ (« - »A)/(t) = (« - ^A) e - ^ ( A + 5 °) t /o(^) (26) 



We can then use an adiabatic approximation to the in- 
tegral (|2"2"j) . which consists in integrating by parts and 
neglecting the second term (which in turn can be used to 
estimate the error in the approximation, as will be shown 
later) : 



,-(A+t(A+* a ))t' 



1 



fo(t')dt' 



A + i(A 



-(A+i(A+5 a ))i' 



X + i(A + S a ) J Q 
The result is then 

K-iA 

<- 2 



fo(t')dt' (27) 



(Ai+i(A + * ))(A 2 + i(A + d«)) 
k — iA 

2gZ~{K-iA){ 1 + i{A + 5 a )) /o(t) 



fo(t) 



(28) 



Finally, this can be substituted back in Eq. (|13bp . which, 
using the explicit definition (IT9l) for /o, yields 



a 



n+ 



2g 2 (n-iA) 



2g^-{n~iA){ 1 + i{A + 8 a )) 



(29) 



The simplest way to handle (f2"9"f is to formally inte- 
grate the right-hand side all the way to the large cav- 
ity roundtrip time (the "quantization time") 2L/c, since 
in that case all the terms in the sum vanish except the 
one with m — n. This may seem unphysical, since, in 
the limit L — ¥ 00, the roundtrip time becomes infinite; 
but in fact what we have on the right-hand-side of (|29|) 
is just something proportional to e ldnt fo(t), and the in- 
tegral of this from to i max will assume its final value 
as soon as i max comfortably exceeds the duration of the 
pulse T, at which point there is no harm in letting ttmax 
formally go to infinity or, in the discrete-mode picture, 
to the roundtrip time 2L/c. Under those conditions, we 
get, for sufficiently large t, 



C n + (too ) — C„+ (0) — 



Cn+(0) 



2 5 2 



2 g 2 («-» A ) M' 22L C fOl 

{K- l A){ 1 + l {A + 8 a )) Mn c Cn+[U) 



2g 2 (n~iA) 



2 k 



2g 2 - (k - iA)(7 + i(A + 8 a )) {ncn/L) 2 + k 2 



(30) 



Again, for a very long pulse (kT 3> 1) it is consistent to replace tick/L in (1301) by A = tiqck/L, in which case, after 
some algebra, the final result can be written as 



C rl +(toc) — 



2 9 2 (k - iA) + (k 2 + A 2 )( 7 + i(A + 8 a )) 
2g\n + iA) - (k 2 + A 2 )( 7 + i(A + 8 a )) 



C n+ (0) 



(31) 



Note that when g = (the pulse does not couple to the atom; for instance, if it has the H polarization when the 
initial state of the atom is zero) one just gets C„+(too) = C„+(0). 
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The result (|3Tj) can also be written with an explicit phase factor pulled out: 

C (t )- K ~ lA y. 2 ff 2 + ( K + * A )(7 + *(A + <U) 

6 " +(<oo) - -^TiA X 2 ff 2-(«_iA)( 7 + i(A + <J a )) C " +(0) (32) 

In this form, the second factor is directly comparable to the reflection coefficient calculated by Waks and Vuckovic 
Q, and used in various limits by many other authors 04^1 • The reason that the prefactor appears here but not in 
the reflection coefficient is discussed in the next section. ^ 



C. Comparison with previous results 

Equation pip clearly shows that when spontaneous 
emission is negligible, the final state of the system con- 
sisting of a single-photon pulse with the "+" polarization 
and the atom in the "+" ground state is merely a phase 
factor times the initial state. This effect has been con- 
sidered before (it is sometimes referred as a "Faraday 
effect" )and quantum logic based on it has also been pro- 
posed pHS! . Any other initial states of the field or of the 
atom, can always written as superpositions of "+" and 
"— " states, from which it follows that the most general 
result one can have is a rotation of the polarization of 
the photon, accompanied by a rotation of the state of 
the atom in the |0),|1) basis. The interesting thing is 
that in general this is an entangling operation. 

More specifically, let 7 = 0. We can write 



Cn+(^oo) 



# C„+(0) 



with 



tan 



(A + S a )(K 2 +A 2 )-2g 2 A 



2g 2 K 



(33) 



(34) 



We also know that C„_ = {Chn — C vn )j\f2 does not 
change with time. We therefore have the transformation 

C hn + C m -e 2 ^ (CUO) + C vn (0)) 

Chn — Cvn — > C/m(0) — C vn (0) (35) 

Adding and subtracting these expressions, we find, in the 
notation of Koshino et al., the following evolution for the 
four possible initial basis states: 

\H,0) -► \H,0) 

\H,l}^ -e^ [i sin <j)\H, 1) + cos <f>\V, 0)] 
|V,0) ->■ -e^ [cos cf)\H, 1) +isin0|V,O)] 
\V, 1) \V, 1) (36) 

When cf> — one has essentially (up to a sign) a SWAP 
gate between the photon and the atom, whereas when 
cf> = 7r/4 one has a VSWAP, as pointed out in To 
have <j> = it is sufficient that all the detunings vanish. 
The condition to have = 7r/4 is, by (j3"4")l . 



(A + 5 a )(n 2 + A 2 ) - 2g 2 (A + k) = 



(37) 



which can always be satisfied, in principle. An impor- 
tant insight of Koshino et al. is that not only is the 



VSWAP a universal gate for quantum computation, but 
by applying it twice one can get an ordinary SWAP as 
well. Hence a VSWAP gate between two photonic qubits 
can be carried out in a completely passive way, without 
any need to directly manipulate the atom, the cavity, or 
the photon detuning, simply by reflecting the first pulse 
twice successively from the cavity, which swaps its state 
for that of the atom; r eflecting the second pulse, which 
carries out a VSWAP between the second photon and 
the atomic qubit; and finally reflecting the first photon 
twice again, so that it acquires the state of the atom after 
the VSWAP. 

The specific results of Koshino et al. can be derived 
easily from (|3T]) by setting 7 = S a — and taking the 
bad cavity limit: k> j,A. Introducing T = 2g 2 / k, one 
has then C^t^) ~ -C„+(0)(r + zA)/(r + zA), and the 
condition for <p = ir/4 is just A = V. One of the key 
results of our analysis is that the VSWAP gate can be 
carried out in a much broader range of regimes, the only 
necessary conditions being that the pulse be very long 
(adiabatic condition) and that spontaneous emission be 
negligible. 

As pointed out in the previous section, many previous 
studies based on the standard input-otput formalism Q 
for the field operators make use of a reflection coefficient 
given by the second factor in Eq. (|32j) . and infer from 
it the phase change of the total atom-field state upon 
reflection. In our formalism, as well as in jL], the total 
phase factor for the state, as written in (|32|) . equals this 
"reflection coefficient" multiplied by a prefactor. We can 
show that if one's specific goal is to compare expectation 
values of the input field operator to expectation values 
of the output field operator, the prefactor cancels, as fol- 
lows: by Eq. (2.6b) of 0|, the output field operator is 



E. 



(+) 



2i ^ 

k 



e a AL 



1/2 



CfeOfce 



-ikL—iQki 



(38) 



The difference between this and the input field operator, 



Eq. flU, is a factor 



in every term in the sum. 



Although one may expect this to be very close to —I for 
all relevant values of k, this approximation is generally 
too drastic, and one needs to use the eigenvalue equation 
([3]) for a better estimate. By writing all the trigonometric 
functions as sums of complex exponentials, it is easy to 
see that ^ is equivalent to 



-2ikL 



-2isinkl - (P/2)e-' lkl 
2ismkl- (i 2 /2)e kl (39) 
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where the approximation f ~ 1 — t 2 /2 has been used 
in the last step. If we then assume kl = 2mr + AZ/c, 
where A, as above, is the incoming field detuning from 
the cavity resonance, and expand on small quantities, we 
can conclude that 



E 



( + ) 



iA 



k — iA 



(+) 



(40) 



In obtaining this result, we used the relationship k = 
ct 2 /4Z between k and the mirror transmission coefficient 
t. Equation (|40| shows that if the state coefficients (|32|) 
are used to calculate the expectation values of the out- 
put field, the prefactor (including the overall minus sign) 
will cancel, and the result will be related to the expec- 
tation values of the input field by only the second factor 
in (1521) . that is, the "reflection coefficient" of the input- 
output formalism. Conversely, note that for reflection 
off of an empty cavity (or one where the atom does not 
couple to the field, g = 0), the input-output formalism 
predicts a phase shift of precisely — (k + iA)/(n — iA) 
(whereas the present formalism yields no change in the 
state vector coefficients); hence, in the input-output for- 
malism, the right-hand sides of Eq. (j3"o]) would all be con- 
sistently multiplied by a phase factor — (n+iA) / (k— iA), 
which, of course, makes no difference physically. Despite 
this agreement, note that, in general, the approach of as- 
cribing to the Schrodinger-picture states the phase shifts 
derived for the Heisenberg-picture field operators works 
only when there is only one photon in the field. When 
dealing with coherent states, for example, there is a big 
difference between the state | — a) and the state — \a). 



IV. NUMERICAL RESULTS 

In this section we show some results of the numerical 
integration of the system (|I3[) in order to examine the 
validity of the adiabatic approximation and to character- 
ize the effect of spontaneous emission on the operation 
of the voWAP gate. Alternatively, one could, in prin- 
ciple, use instead the exact (in the continuous limit) so- 
lution (1221) . which, for the Gaussian pulse we will con- 
sider, can be evaluated in terms of error functions of 
complex arguments; but, after that, the algebra becomes 
cumbersome, whereas the direct numerical integration of 
Eqs. (fT3")) poses no particular challenges. It is, however, 
necessary to realize that, even if the initial pulse has a 
very small bandwidth, once it is inside the cavity this 
changes, and one may need to consider a very large num- 
ber of "modes of the universe," of the order of kL/c or 
TL/c, in order to accurately determine the asymptotic 
behavior investigated below. This point is elaborated in 
more detail towards the end of this Section. 

We take the spectrum of the initial pulse to be given 
by the coefficients (fT2")l . out of which we can build the 
initial coefficients C„+(0), which will formally be equal 
to C n /^2. Then C n _(0) = ±C n /V2, depending on 
whether we take this to be a vertically or horizontally 




FIG. 2: Fidelity F versus coupling constant g for an initial 
Gaussian pulse of duration T = 100/k;. For every value of 
g/n, A/k is chosen so as to satisfy the condition Eq. (I37[) . 
with 8 a — 0. 



polarized pulse. In any case, the C n ~ (0) are constant, so 
for a measure of the performance of the gate it suffices 
to consider the C n +(0). In particular, we shall define 



Fe ** = E„Q + (o)c n+ (oo) 



£„C* + (0)C„ + (0) 



= -2^C; i+ (0)C7„+(co) 

n 

(41) 

where F is a measure of the "fidelity" of the final 
wavepacket to the initial one, and $ measures the phase 
shift of the state in the Schrodinger picture. Comparing 
to Eq. ([33"]) we see that if the adiabatic approximation 
were exact, one would have F = 1 and $ = <f> as given 
by Eq. (f34|) ; the deviation from these values can be used 
to characterize the gate error. 

To begin with, we show in Fig. 2 that the gate works 
well in both the bad cavity and good cavity limits, pro- 
vided only that the pulse is sufficiently long. This re- 
quires both T > 1/k and T > 1/T (with T = 2g 2 /n). 
For the parameters of the figure, the first inequality is 
always satisfied, whereas the second one requires g/n ^> 
0.07. Note that this includes both bad cavity (g < k) 
and good cavity (g > k) values. 

The way the gate fidelity improves as the pulse dura- 
tion increases is illustrated in Figure 3, for V = 0.5k (i.e., 
g = k/2) and variable T. The best fit slope in the log- log 
plot is —1.99, which suggests that the fidelity improves 
quadratically as T increases. 

This quadratic dependence may be a little surprising, 
since a glance at Eq. ([27| suggests that the terms ne- 
glected in the adiabatic approximation only decay as 
1/T. This is true in general, but a closer examina- 
tion reveals that for a long pulse with a symmetric fre- 
quency spectrum, the more favorable 1/T 2 can be ap- 
proximately realized. To see this, note that an adiabatic 
approximation to the second term in (1271) yields essen- 
tially e- i ( A+5 ») t / (i)/[A + i{A + 5 a )] 2 , where by Eq. (USD 
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FIG. 3: "Infidelity" 1— F versus pulse duration T for g — k/2, 
and A and S a as in Fig. 2. 



FIG. 4: Difference between the phase factor $ and 7r/4 versus 
pulse duration T for the same parameters as in Fig. 3. The 
best straight-line fit has slope —2.17. 



/o goes as 

fo(t) = ie« A + 5 ^ J2 <Cn+(0) (A + 5 a - S n ) e i<5 "' 

n 

(42) 

When this is substituted in Eq. (|13bj) and the integration 
over the quantization time is performed, one obtains a 
correction to C„+ which goes as 



(nor / L) 2 + k 2 

n — no 
(nor / L) 2 + k* 



C n+ (0) 



(43) 



(making use of the fact that 5 n = c(n — no)ir/L + A + 
<5 a , as indicated below Eq. pi))). For a pulse such as 
(fT2")l . for which |C„| 2 is an even function of n — tiq, the 
sum ~^2 n <5C n +(oo)C* + (0) will vanish, which means that 
the next-order contribution to the infidelity will come 

from /o ~ 1/T 2 , in agreement with the numerical results. 
Figure 4 shows that this also applies to the gate phase 
error, that is to say, the d ifference between $ and the 
value 7r/4 needed for the VSWAP gate also appears to 
decrease at least quadratically, although here the results 
of the numerical integration are not so clean. 

The theoretical possibility to carry out a determinis- 
tic gate between single photons with an error scaling as 
the inverse square of the gate duration is remarkable in 
view of the results derived by one of us in [l(| for quan- 
tum logic with atomic qubits. There it was shown that 
if the gates were mediated by a quantized field in a co- 
herent state, with an average number of photons n, the 
gate error would scale as 1/n. Note that this constraint 
would apply to the original Duan-Kimble gate 2] , which 
requires the manipulation of the atom by a coherent field 
in between single-photon pulse reflections; in contrast, 
the Koshino-Ishizaka-Nakamura gate appears to be able 
to work with negligible error at the single-photon level. 

The 1/T 2 scaling derived here is also better than the 
1/T scaling postulated to hold, alternatively, in [l(| for 



quantum logic with material qubits and quasistatic con- 
trol fields. Time and again over the past few years one or 
the other of these constraints has been shown to hold for 
all kinds of systems of experimental interest [111 ], and in 
many cases it has been shown to follow from very funda- 
mental considerations, such as spontaneous emission [l2j], 
or conservation laws [13l. Il4j. The present result, how- 
ever, appears to suggest that quantum logic should be 
possible with a much smaller energy cost than had pre- 
viously been thought possible (although one cannot yet 
rule out the possibility of some hidden energy cost in, 
for instance, the electronic devices that might be used to 
route the photons towards or away from the cavities). In 
any case, this is a question that certainly requires further 
investigation. 

Besides the breakdown of the adiabatic approximation, 
an important factor that may degrade the performance 
of the gate is spontaneous emission. Figures 5 and 6 com- 
pare the fidelities and phase shifts obtained numerically 
for zero and nonzero 7, and show that the damage done 
by 7 is approximately linear for the fidelity and quadratic 
for the phase angle. These scalings agree with the result 
of a simple expansion in powers of 7 of the adiabatic 
approximation Eq. (pTTj) . 

We note that our treatment of spontaneous emission is 
formally equivalent to considering only the first term in 
a "stochastic wavefunction" unraveling of the underlying 
master equation [l5l - tl8j . Up to an overall normaliza- 
tion factor, this just gives the evolution of the system 
conditioned on the probability of no actual spontaneous 
emission event taking place. The probability to, in fact, 
lose the photon because of spontaneous emission can be 
estimated from the adiabatic approximation as 



P loss = 2 7 / \C 2 {t)\ 2 dt 
h 



(2g 2 -A(A + S a )) 2 + K 2 (A + S a ) 2 



(44) 
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FIG. 5: Difference between the fidelity obtained numerically 
when 7 = and when 7 is nonzero. The pulse duration 
T — W0/k, and other parameters are as in Figs. 3 and 4. 
The best-fit slope is 0.98. 
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FIG. 6: Difference between the phase factor $ obtained nu- 
merically when 7 = and when 7 is nonzero. Parameters as 
in Fig. 5. The best-fit slope is f .94. 



Making use of the condition (|37| which the detunings 
have to satisfy in order to have a VSWAP gate, this 
simplifies substantially, to 



K 



(45) 



which suggests that the optimal detuning arrangement, 
to minimize Pi OS s, is to have A = (pulse resonant with 
the cavity) and <5 a = 2g 2 / k (atom detuned from both the 
cavity and the pulse). In that case, the probability that 
the photon be lost through spontaneous emission is just 
the inverse of the "Purcell factor," F — 2g 2 /nj. 

Finally, it is also possible to use the numerical results 
for C n + if) to visualize the evolution of the pulse in space, 
as well as in the frequency domain. Figure 7 shows the 
pulse intensity I- m (t,z) — (E^'(t, z)E^\t, z)) at three 
instants, before, during and after the interaction with the 



FIG. 7: Pulse intensity profiles, at three different times, for 
the case of a vertically-polarized pulse incident on the cav- 
ity under VSWAP conditions. The reflected pulse has both 
horizontal- and vertically-polarized components. The posi- 
tion z has been scaled by the distance that the pulse travels 
in a cavity lifetime, that is, z s — zk/c. The dotted verti- 
cal line represents the position of the cavity and the view is 
"unfolded" so that z a > corresponds to a reflected pulse 
(see text for details). The parameters for the calculation are 
T = 0.1k, T = 100/k. 



cavity; the initial field-atom state is taken to be | V, 0) , 
and the detuning is taken to satisfy the condition for the 
%/SWAP, in which case Eq. (|5o| predicts, and the nu- 
merical calculations show, that the initial pulse intensity 
is reduced by a factor of 2 and a horizontally polarized 
pulse with the same reduced intensity is generated by the 
interaction. Note that, by Eq. and the discussion fol- 
lowing it regarding the phase factors e 2lkL , the outgoing 
field intensity I ut(t, z) is formally identical to I\ n {t, —z), 
so the part of Figure 7 corresponding to z > can be 
regarded as an "unfolded" (reflected around the z = 
plane) view of the evolution of the outgoing pulse in the 
one-sided cavity setup of Figure 1. This also means that 
the cavity itself is invisible in the figure; the vertical line 
at z = merely represents its input (which is also its 
output) mirror. 



The need for the spectral coefficients C„+ to reproduce 
the sharp discontinuity exhibited in Fig. 7 by the analyt- 
ical functions Ih(t, z) and I v (t, z) explains why it is neces- 
sary to keep a large number of coefficients in the calcula- 
tion for accuracy. Figure 8 shows that the spectra them- 
selves are clearly broadened somewhat when the pulse is 
in the cavity; in particular, a close scrutiny shows that 
there are long, but relatively low magnitude, tails that 
cover the whole cavity bandwidth. Figure 8 also shows 
that the VSWAP condition, here satisfied by S a = and 
A ~ r, results in some frequency pulling of the spectrum 
towards the cavity line center (at n = 0, or n—n = —no), 
while the pulse is in the cavity. 
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FIG. 8: Frequency spectra of the different polarizations for 
the situation depicted in Fig. 7. 

V. CONCLUSIONS 

In this paper, we have extended the analysis of [l[ to 
show that the kind of all-passive, cavity-mediated, single 
photon quantum logic proposed by Koshino, Ishizaka and 
Nakamura actually can work in the good cavity limit as 
well as in the bad cavity limit originally considered in [1| . 



We have derived formally exact analytical expressions for 
the excitation probability, and approximate results valid 
for sufficiently slow pulses. We have shown that the gate 
error probability in the absence of spontaneous emission 
scales as 1/T 2 , where T is the pulse duration, and given 
expressions for the photon-loss probability due to spon- 
taneous emission valid also in the adiabatic limit. Our 
formalism also has allowed us to clarify the relationship 
between the phase factor for the reflected pulse derived 
by Koshino et al. and previous results for two-level sys- 
tems found in the literature. 

We believe that the possibility of deterministic quan- 
tum logical gates at the single-photon level, without the 
need for complicated manipulations, opened up by the 
work in [![, deserves a great deal of attention. Many 
groups are currently working on possible schemes for all- 
optical quantum computers [l9| that suffer from enor- 
mous overheads because of the probabilistic nature of 
available photon-photon logical gates. Even if an all- 
optical quantum computer turned out not to be the best 
choice in the long run, deterministic gates for single pho- 
tons could still prove invaluable for quantum communica- 
tion, to perform, for instance, entanglement distillation 
or state purification on photons carrying quantum infor- 
mation across nodes of a quantum network. Further stud- 
ies are clearly necessary to identify suitable systems, and 
to characterize the performance of these types of gates in 
realistic setups with currently- available technology. 

Finally, as we have indicated in the previous section, 
the possibility of these gates also raises questions of fun- 
damental interest, such as the true minimum energy re- 
quirements for quantum computation, which we intend 
to pursue in the near future. 
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